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Abstract 



■ We consider families of dynamics that can be described in terms of Perron-Frobenius 

operators with exponential mixing properties. For piecewise expanding interval maps 
we rigorously prove continuity properties of the drift J(A) and of the diffusion coefficient 
• D{X) under parameter variation. Our main result is that D{X) has a modulus of continuity 

' of order £9(|(5A| • | log |5A|)^), i.e. D{X) is Lipschitz continuous up to quadratic logarithmic 

corrections. For a special class of piecewise linear maps we provide more precise estimates 
at specific parameter values. Our analytical findings are verified numerically for the latter 
class of maps by using exact formulas for the transport coefficients. We numerically observe 
strong local variations of all continuity properties. 

> 

1 Introduction 

(N ■ In simple deterministic dynamical systems physical quantities like transport coefficients can be 

fractal functions of control parameters. This finding was first reported for a one-dimensional 
C3 I piecewise linear map lifted periodically onto the whole real line, for which the diffusion coefficient 

was computed by using Markov partitions and topological transition matrices |26l[?7l[^ . A 
generalization of this result was obtained for a map with both drift and diffusion by deriving 
exact analytical solutions for the transport coefficients [151 15] . Further maps modeling chemical 
reaction-diffusion 14J and anomalous diffusion [33] yielded also fractal transport coefficients. 
}_( ' Recent work aimed at physically more realistic models like (Hamiltonian) particle billiards, 

for which computer simulations yielded transport coefficients that are non-monotonic under 
parameter variation [31] . Ref. [32j contains a summary of this line of research. 

These results asked for a more detailed characterization of the "fractality" of transport 
coefficients. A first attempt in this direction was reported by Klages and Klaufi f30], who 
used standard techniques from the theory of fractal dimensions for characterizing the drift and 
diffusion coefficients of the map studied in [15j . They numerically computed a non- integer 
box counting dimension for these curves which varied with the parameter interval, leading to 
the notion of a "fractal fractal dimension" . These results were questioned by Koza [35] , who 
computed the oscillation of these graphs at specific Markov partition parameter values. His 
work suggested a dimensionality of one by conjecturing that there exist non-trivial logarithmic 
corrections to the usual power law behaviour in the oscillation. 

This research reveals the need to study the parameter dependence of transport coefficients 
in a rigorous mathematical setting, which can be formulated as follows: Given a parametrized 
family of chaotic dynamical systems T\ : / ^ / on an interval / with unique invariant physical 
measures fix together with a family of sufficiently regular observables ipx : I R one has, under 
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suitable mixing assumptions on the systems {T\, ^x), a law of large numbers and a central limit 
theorem for the partial sum processes S\^n{x) = ^22=0 V-'a(T'a^), namely 

lim n^^Sxn — JW '■— / dn\{x) for /^(-a.e. x 

n^oo ' Jj 

and 

where D{X) :— lim„^oo ^ Jj (X]fe=o('/'A(^A^) ~ ^Wij diix{x). For suitable choices of the 
observables ipx, the process Sx,n is just the deterministic random walk generated by a lift of 
the map Tx to the real axis, and J(A) and D{X) are the drift and diffusion coefficient of this 
random walk respectively. 

There are a few rigorous results in the literature describing the dependence of /Jx and of 
quantities like J(A) for various classes of systems. Without going into the details they can be 
summarized as follows: If the maps Tx and the observables ipx depend smoothly on A and if 
the topological conjugacy class of Tx is not changed when A is varied, then ^x (and hence J{X)) 
depends differentiably on A [ill51ITUllTOll20ll39HD] . If the topological class changes, quantities like 
J(A) may behave less regular and have a modulus of continuity not better than \SX ■ log\SX\\, 
even for very simple maps Tx like symmetric tent maps [3 . On the other hand, this modulus of 
continuity is the rule for systems whose Perron-Frobenius operator (acting on a suitable space 
of "regular" densities) has a spectral gap [22l[24]. 

The goal of this paper is to explicitly relate these mathematical results to transport coefh- 
cients. We do so by rigorously proving continuity properties of J(A) and D{X) under parameter 
variation for certain classes of deterministic maps. In Section [2] we give a general estimate 
for families of dynamics (deterministic or not), which can be described in terms of Perron- 
Frobenius operators with exponential mixing properties. The applicability of these general 
results to piecewise expanding interval maps and in particular to the class of piecewise lin- 
ear maps discussed in t26,.27,J9,,_15,_32j is checked in Section [3l The main result is that D{X) 
has a modulus of continuity of order Od^Aj • | log |(5A|)^), i.e. D{X) is Lipschitz continuous up to 
quadratic logarithmic corrections. In Section 2] we summarize the general results for transport 
coefficients in the special case of piecewise linear maps and provide more precise estimates for 
special parameters. Our analytical findings are verified by numerical computations in Section[5l 
for which we use exact analytical formulas of the transport coefficients [15]. Particularly, we 
numerically analyze local variations of these properties. Our work corrects and amends previous 
results reported in [501155] . 



2 The general setting 

Let / be a compact interval, m normalized Lebesgue measure on /, the space of Lebesgue- 
integrable functions from / to R, and BV C £^ the space of L^-equivalence classes of functions 
of bounded variation. We use the following simplified notation for the two corresponding norms: 

l/li := J I/I dm, 11/11 :=Var(/) (1) 

where 

Var(/) := sup f^' dm : ^ e C\R,R),\^\oo < 1^ (2) 

is the variation of / as a function from M ^ M (i.e. extended by / = on M \ /). If / 
is differentiable as a function from K ^ K integration by parts shows easily that Var(/) — 
J l/'l dm. Var is obviously a semi-norm, and as |/|i < |/|oo < i Var(/), it is actually a norm. 
This and more details on functions of bounded variation can be found in [25l section 2.3]. The 
monograph ^ is a comprehensive reference for most of the background material needed in this 
section. 
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We consider a family T of nonsingular maps T : I ^ I. Nonsingular means that the 
Perron- Frobenius operator Pt '■ L\-^ L\-^ is well defined, i.e. 

j Prf ■ gdm ^ I f-goTdm (/ e L^^ge L^) . (3) 

By definition, |Pt|i = 1 for all T E T, and we assume 
Hypothesis 1 Ci := sup{||P-P|| : T e T, ?i e N} < oo. 

Our main assumption is that the maps in T are uniformly exponentially mixing in the following 
sense: 

Hypothesis 2 Each T G T has a unique invariant probability density hx G BV (so Pxhx — 
hx), and there are constants 7 G (0, 1) and C'2 > such that, for all T G 7", 

I^t/Ii < (^27" Var(/) for all f e BV with J / dm = and for allneN . (4) 

Observe the following consequences of Hypothesis [1] and [2] 

iPrf - hrli < C27"(Var(/) + 2Ci) for all probability densities/ G BV (5) 

and 

Var(/iT) < 2Ci {T € T) . (6) 

Indeed, |P^/-/it|i = |P?(/-/it)|i < C27"(Var(/) + Var(/iT)) ^ as n ^ 00 by Hypothesis^] 
for each probability density / G BV, and Var(Ppi) < Ci Var(l) = 2Ci by Hypothesis[TJ Hence 
^ follows from the definition ^ of Var(.), and then ^ is an immediate consequence. 

Since it is our goal to investigate the dependence of various dynamical quantities as functions 
of T G T, we need to introduce a distance on T. At this stage the following one, which was 
already considered in [22j . is most apropriate. It measures the distance between two maps Ti 
and T2 from T in terms of a suitable norm of P^i — • 

III^Ti - PtJII sup{|Pti./ - PtJ\i ■ f ^ BV, \\f\\ < 1} . (7) 

This distance can be controlled in terms of a more "hands-on" distance between the graphs of 
the maps: 

d{Ti, T2) := inf{e > : 31^ C / and 3 a diffeomorphism a : I ^ I s.th. 

m{I \ h) < e. Til/, = T2 o a\i^ , and (8) 
Vx G /e : \aix) ~x\<e, \1/(t'{x) - 1| < e} . 

Namely (see [211 Lemma 13]), 

III^Ti -PtJII < 12-d(ri,T2) . (9) 
Now, as a warm-up exercise, we can prove the following estimate: for fc > let 

4:(0,cx3)~.(0,oo), 4(«) :=«•(! + I logul)^ (10) 



Lemma 1 There exist constants K'l , Ki > such that 

\hT,-hT,\i<K[-h{\\\PT,-PTA\)<Ki-li{d{Ti,T2)) (ri,T2Gr) (11) 
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Proof: Let fj :— \\\Pti — Pt2\\\^ assume without loss of generality that r) < 1, and fix e N. 
For / e BV, 



\PtJ - PtJ\i < E K^'^'iPT. - PT.)Py\i < E - Pt.)PtJ\i 

k=0 k=0 

< E III^T, - PrMP^fW < E ^C'lll/ll < CiiVr?||/|| 



(12) 



fc=0 fc=0 

where we used Hypothesis [TJ Hence, 

\hT. - /^T.ll < - P^lll + K(l - hT,)\l + |P^(1 - /itJIi 

< 2CiNr~i + 2 ■ C2l^ (2 + 2Ci) 
where we used Q and ®. With N = \^^], this is (HH). q.e.d. 

Remark 1 Even if T is a family of piecewise linear maps and if Ti has the Markov property, this 
estimate can generally not be improved. Examples for this fact within the family of symmetric 
mixing tent maps are provided in [3l[37]. 

Suppose now that to each T £ T there is associated an "observable" ipr : / ^ M. We make 
the following assumptions: 

Hypothesis 3 C3 := sup{Var(i/iT) : T e T} < 00 

Hypothesis 4 There is C4 > such that \iPt, - Vt^ |i < C4d(Ti, Tz) for all Ti, T2 G T. 



Denote 



J(T) := J iPtHt dm 
Then we have immediately from ^ and Lemma [T] 



(13) 



Corollary 1 There is some K2 > such that, for all Ti,T2 G T , 

\JiTi)-J{T2)\<K2-ii{diTi,T2)) (14) 
J(r) is the "drift" of the partial sum process 



ST,n E ° = nJ{T) + Y,^ToT'' 



k=Q k=Q 



under the invariant measure hxm, where V't = "tpT — J{T). Observe that 

Var(^T)<2C3, l^n - VtJ < 2C4 d(Ti, Ts) for all T, Ti , T2 G T . (15) 
In view of Hypothesis [5] we can also define the "diffusion coefficient'll] of this process: 



D{T) := lim 



— ( 'tp^ hx dm + / 'ij^T • V'T o /iT t^'TT- 

— ( -tp^ hx dm + / PxitpThT) dm 

J n=l"' 



(16) 



^This is the convention in the physics literature. In the mathematics literature one would rather call 2D(T) 
the diffusion coefficient. 
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Even more, we have the central Hmit theorem 

C{n-^ (Sr^ri - nJ(T)) ^ JV{0,2D{T)) as n ~* oo , (17) 

see e.g. [2TJ[T8l[38]. Among physicists (fT6|) is known as the Taylor-Green-Kubo formula for 
diffusion [32]. For the dependence of D{T) on T wc prove: 



Proposition 1 There is some > such that, for all Ti,T2 G T, 

\D{T^)~D{T2)\<K^-t2{d{T^,T2)) 
Proof: Observe first that, for all ip,h E BV, 

Var(V'/i) < Var(V')|/i|oo + Var(/i)|V'|oo < Var(V') Var(/i) 



(18) 



(19) 



Indeed, for differentiable u we have Var(u) = / \u'\ dm, so for differentiable h and (fT9|) follows 
from the product rule of diferentiation. General h and are then approximated using moUifiers. 
It follows that, in view of 



P^lipxhx) V'T dm 



< \P¥i^^ThT)\l ■ I^tIoo < C2C37"Var(V;T/iT) 



(20) 



Let 77 = ||Pti - -Pts II as before, denote 77 := d{Ti,T2) (so that 7) < 12tj), and fix iV = \^^^■ 



For all Ti,T2 G T, eq. 1^ implies 



Pti (VTi /iTi ) VTi dm- P^^ (ipT^ ) i/Ts 



16C1C2CI 
1-7 



For < n < N we use a different estimate. We decompose 

Pt, i^T, hT, ) iPt, dm- [ P^^ {tJjt, Ht, ) Vt. dm = A^' + + + 



where 



|A?| := 



Pti (V^Ti /iTi ) (VTi - ■4'T2 ) 



<ivar(P^^(V'T,/iTj)|^T, -^tJi 



(21) 



(22) 



(23) 



and 



|A^| 



/ {P¥,-P¥J{4'nhT,)i^T,dm 
<ACiCl\\\P^^ - P^J\ < 4(CiC3)'n77 
where the last inequality follows from eq. Next, 



< |||P^^, -P^,|||Var(V'T,/iTj|V'Tjc 



(24) 



|A^| := 



Pt2 {ii'Ti - ^T2)fT'Ti)4'T2 dm 

<ACiC3.Ciri 



< \lpT^ - V^tJi |/lTi|oo IV'Talc 



(25) 



and 



Pt2 (^T2 (/iTi - )) V'Ts 



< |/lTi - ^Tall l^/Taloo I^/TjIc 



(26) 



<CfXi4(?7) 
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^From (1231) - (EH) we see that 

\A-\ + \A^\ + \A-\ + \A-\<K{nr^ + iM) (27) 
for some constant K > 0. Hence, in view of pi]) and the choice of N, 

^ n— 1 

AC.r^r'^ (28) 

< V + K{N'v + (2iV - l)£i(77)) 

<i^3-^2(r/) 

for a suitable constant K^. q.e.d. 



Remark 2 Quite often slightly stronger forms of Hypotheses HI and HI are satisfied, where the 
mixing assumption ^ is replaced by 

Var(P^/) < C27" Var(/) for aU f £ BV with J / dm = and n e N (29) 

and the assumption on the T-dependence of '0T is strengthened to 

Var(^Ti - V'tJ < Cid{Ti , T2 ) for aU Ti , G T. (30) 

An inspection of the above estimates shows that | A"| < ACiO^C^Ci'j^r] and | | < 4CiC2C3(l— 
7)~^77 if is assumed. If additionally (|30p is assumed, then IA3 | can be estimated as follows: 
Let a := Jj{iJjti — ipT2)hTidm. Then 



J ((VTi - V'T2)^Ti - a/iT2)V'T2 < ^2(^3 7" Var (^{ipTi - '0T2)^Ti - a/iT2) 



(31) 

Hence, + l^al = (^(^) uniformly in N and = 0(A^^) = 0(^1(77)). But 

we see no way, in general, to bound the AJ-terms in a similar way. However, for particular 
families of maps (which are all topologically conjugate), we will see in subsection 14.21 that 
A J — for all n and that the estimate for IA2 | can be made more precise. 

3 Checking Hypothesis [1] and [2] 

3.1 General piecewise expanding maps 

In this subsection we show how the general Hypothesis H] and 121 can be verified in the more 
particular setting when T is a parametrized family of piecewise twice continuously differentiable 
and expanding interval maps. So, from now on, we look at the following setting: 

A C M'' is a compact parameter space, T = {Tx : A G A}, and (Tl) 
there is some L>0 such that d^Tx^^Tx.^) < L\Xi - A2I for aU Ai, A2 G A. (T2) 

We start with an abstract result which reduces Hypothesis lH essentially to a uniform Lasota- 
Yorke type inequality. 

Lemma 2 Assume (|Tip and (|T2p . Then Hypothesisl^ andl^ are valid if the transformations 
T £ T are mixing and satisfy a uniform Lasota-Yorke type inequality.' there are constants 
C'5, Cg > and a G (0, 1) such that 

Var(PJ,7) < Cga" Var(/) + Cel/li for all T e T,n e e BV. (LY) 
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Proof: As |Pt|i = 1 and |/|i < ^H/H, it is straightforward to check that Hypothesis [T] holds 
with Ci = C5 + iCg. 

We turn to Hypothesis [H Note first that, because of (jTip and (|T2I) . it suffices to show that 
for each A G A there are (5(A) > 0, C2(A) > and 7(A) G (0, 1) such that ([4]) holds with these 
constants for all Txi with |Ai — A| < (5(A). But this is guaranteed by Corollary 2 and Remark 
Ic in [21]. q.e.d. 



Our next task is to give sufficient conditions for (jLYp . To this end we specialize further and 
assume from now on that our maps are piecewise expanding (PE) maps in the following sense: 

For each A G A there is a finite partition , . . . , I^^ ) of / into sub- 
intervals such that all T\\jj are monotone, C^, and k\ :— inf \T^\ > 2 . ^ 

Already in [36] it was proved that each individual (PE)-map (even if 1 < ka < 2) satisfies 
(|LYp with constants C5,Ce,a depending on the map. For parametrized families of maps one 
can generally find uniform constants, but there are counterexamples where this is not possible 
|22[ [6l[7]. Under the above assumption ini\^\K\ > 2 one can, however, give simple sufficient 
conditions ensuring the uniform LY-inequality. The proof in (see also [25l Proposition 2.1]) 
shows 

Var(PT./) < — Var(/) + {Dx+Ex)\fh (32) 



where 

Dx = sup 



K\ min," 1/^ ' 



(33) 



From this (jLYp follows with a ~ = 1, and Cg = sup;^ 1^^^^^^ + -^^) provided this 

supremum is finite. 

In some cases of interest the i?A, A G A, are not bounded because there are arbitrarily short 
monotonicity intervals. In such situations, ad hoc arguments are needed. We give an example 
in the next section. 

3.2 Piecewise linear modulo 1 maps 

We now look at a particular model dealt with in [TWIMTB1IT7] from a mathematical perspective 
and in [ ^ [?7 l l ^ n ^ l32j from a physics point of view. Let / = [-i, i], A = [ao,ai] x [-i, \] 
for some constants 2 < oq < oi, and for A = (a, 6) G A consider 

Tx{x)^ax + b mod(Z~i). (34) 

Hofbauer [16] showed that these maps have always a unique invariant probability densitjH, but 
although these maps received further attention also in the mathematical literature [T7lll2lll3j , it 
is not so easy to draw Hypothesis [2| from these sources. Therefore we will take up a rather direct 
computation made in [23] to prove, without having to rely on the compactness assumption (jTip . 
the following lemma. 



Lemma 3 Let A — (a, h) G A. Then 

YaiiPT^ f)<^ Var(/) + 2 / / dm 



for all f eBV andne N. (35) 



(This implies immediately Hypothesis Ui with Ci = 1 + and Hypothesis \^ as well as its 

strengthening ([29]) with 7 = ^ and C2 = C2 = l.J 



■^Indeed, Hofbauer shows this for the maximal measure of such maps, but since these maps have constant 
slope, the maximal measures are just the absolutely continuous ones. For numerical results on the probability 
densities associated with these measures and how they change under parameter variation see 1271 . 
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Proof: Denote by T the family of all -functions : M ^ K with \(p\oo < 1- In [131 eq. (11)] 
a number T{ip, s) > is defined for each pair of if d J- and s G M. In view of 23j eq. (13) and 
(14)] it suffices to show that for each ip £ T there is some s € [—2, 2] such that T{(p, s) < 
We will show that this is the case for s = s,^ fi'^) ~ (p{—^)- 

Let p = - 6] - 1 and q = + 61-1. Then -p - i < Tx{~^) < ~P + \, 

g— i<Tx(i)<q + i, and Tx has monotonicity intervals (A^, 3^), fc = — p, . . . ,q, where 

A-p = -^, Ak^a-\k-^-b) {k = -p+l,...,q) 

Bq^\, Bk^a-'{k+^-b) (fc = -p,...,g-l) 

In order to estimate r((^, s) in [531 eq. (11)] one has to evaluate certain terms Uk ■— ^{Ak) — 
h9k{Ak) - sAk and Vk := '0(^fc) - hgk{Ak) - sBk- In our case, gk{Ak) = a~\ as gk is the 
inverse of the derivative of the fc-th monotone branch. The quantity ik is an abbreviation for 
ip{TxAk), where TxAk denotes a limit from the right. Therefore, using the formula on the 
bottom of [531 p. 1779], we obtain 

aUk^a- (ipiAk) - ikgkiAk) - s^Ak) 

= - fi-^ + b + p) + (fc +p)((^(i) - (^(-1)) _ . (fc - 1 - 6) 

a 1 (36) 

= -'p{-- + b + p) + s^-{p+- + b) iffc>-p, 

a U-p = - (p{-- + b + p) + s^- - , 
and similarly, 

aVk =a • {ij{Ak) - ikgkiAk) - s^Bk) 

^-v{-^ + b + p) + s^-{p-'^ + b) iik<q, (37) 

aVg=- v{-^+b + p) + ■ {p + q - ^) . 

It follows that \Vk — Uk'\ < a^^|s;^| for all k,k' e {—p, . . . ,q}. Hence, by [23l eqs. (11) and 
(12)], 

r{ip, s^) < fl-i + -a-^\s^\ < - . (38) 

q.e.d. 

Next we check assumption (jT2l) on the Lipschitz dependence of the maps on the parameters, 
so we estimate diTa.b^Ta' ^b')- For the proof we extend the maps to the whole real line (keeping 
the same names) by applying definition (I34p to all a; G M. 

Suppose a' < a and denote by Ak and A'j^, k = —p + 1, . . . ,q, the discontinuity points of 
the two maps as introduced in the proof of Lemma [S] Consider the linear map L : R ^ K, 
L{x) = (ax + b — b')/a' and observe that L{Ak) — A'j^ {~p < k < q) and a'L{x) + b' — ax + b 
for all x S M. Let [u, v] := / n L^^{I) and /o '■— [u + S,v — S] for some arbitrarily small 6 > 0. 
Define ct : / / by <y{x) = L{x) if a; G /q and extend cr to a diffeomorphism of /. Then 

• m{I \ /o) < (1 - ^) + \b' - b\/a + 26 < flg i(|a' - a\ + \b' - b\) + 26 

• \(j{x) ~x\^ \L{x) - a;| < - 1| + |6 - b'\/a' < \i|a - a'\ + \b - b'\) for all x £ lo 

• \1/<t'{x) - 1| = - 1| < flo V - a\ for aU x e Iq 

• Ta'(a(a;))=r„(a;)forallxe/o. 

So d{Ta,Ta') < ^(|a — a'\ + |6 — b'\) as 5 > could be chosen arbitrarily small, 

'^var(.) in |23l is the same as Var(.) here. This is different from the use of var(.) in |25l . 
Following this reference precisely, the reader will notice that instead of the factor ^ one gets the factor 
^ + ^Vg. But Vg = for piecewise linear transformations as noticed at the bottom of '23' p. 1779]. 
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4 Transport coefficients 



We apply the results of the previous sections to determine transport coefficients of the deter- 
ministic random walks generated by the maps Tx — Ta^ from subsection 13.21 The random 
walks in question are Sx^n = Sfc=o ''/'^ ° -^a with 

ijx{x) ^ {a - l)x + b . (39) 

It is an easy exercise to see that Hypothesis [3] and |3] as well as their strengthening ([50]) are 
satisfied: Var(^/>a,fc) < 2(a-l + |6|) < 2(ai + l) =: C3 and \ipa,b-'4'a' ,b'\i < 5 Var(V'a,b-''/'a'.b') < 
\a-a'\ + \b-b'\. 

For later use we note that the maps Ta,b and Ta,^b are conjugate in the sense that Ta,fc(— a;) = 
^Ta.-b{x), in particular ha.-b{—x) is also an invariant density for Ta.b and, by uniqueness, 

ha,b{,x) = ha-b{-x). 

We first note the following explicit form of the drift: 



J(A) = J(a, b) :=: j ipa.bha.b dm = b + {a - 1) J xhafi{x) dx . (40) 
As noted above, hafiix) — hafii—x). Hence J(a,0) = oH 

4.1 Upper bounds for the modulus of continuity of the drift and the 
diffusion coefficient 

Now we apply Corollary [1] and Proposition [1] to our setting. 

Proposition 2 For the family of maps (Tx : A € A) defined above, there are constants K^, K4 > 
such that the drift J(A) := J{Tx) and the diffusion coefficient -D(A) := D[Tx), A = {a,b), 
satisfy 

|J(A)-J(A')|<i^3-|A-A'|-(l+|log|A-A'||) (A,A'eA) (41) 
|i?(A)-i?(A')| <i^4-|A-A'|-(l+|log|A-A'||)' (A,A'eA). (42) 

Corollary 2 a) The graph of D : A ^ W has box- and Hausdorff- dimension 2. 

b) For each 6 G M, the graph of Dh : [ao,ai] — > M, Di,{a) = D{a,h), has box- and Hausdorff- 
dimension 1. 

c) For each a > 2, the graph of Da '■ [~^, ^] ~> K,. Da{b) = D{a,b), has box- and Hausdorff- 
dimension 1. 

Proof: Denote by dims and dim/f the box and Hausdorff dimension, respectively. Obviously, 
2 < dimij(graph(D)) < dimB(graph(£')). So it remains to show that dimB(graph(iI')) < 2. To 
this end subdivide the rectangle A into little squares of equal size « . For each such square 
Q we have 

max{i:)(A) : A G Q}-min{D{\) : A G Q} < K -'^ [l + log Nf . (43) 

Hbiicc 

dim.(graph(i^)) < limsup ^^^i^^^^i^ ^^^^^ ^ ^ 
The two other claims are proved in the same way. q.e.d. 



^In 1271 it was conjectured, based on analyzing these deterministic random walks in terms of Markov par- 
titions, that for 6 = and a > 2 the maps Tx exhibit a central limit theorem and that a diffusion coefficient 
exists, which is confirmed by II17I I. 
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Remark 3 Corollary 2 has already been conjectured by Koza [5S]. His conjecture was based 
on calculating the pointwise Minkowski-Bouligand dimension for algebraic Markov partition 
parameter values of this family of maps by using the exact solutions for drift and diffusion 
coefficient given in [15]. This led him to conclude that the oscillation [41] of D{X) is linear in 
the size of the subinterval multiplied with a logarithmic term, cp. (20) of [3S] with above. 
The exponent of this logarithmic correction was found to be either one or two depending on 
the type of Markov partition. 



4.2 A closer look at maps with integer slope 

We finish this section with a closer look at the functions Da{h) when a is an integer larger 
than 2. In this case T = Ta.b can be seen as an a-fold covering linear circle map, so it leaves 
Lebesgue measure invariant. Therefore = 1 for all such T and the estimates from the proof 
of Proposition [1] simplify drastically: Fix a S {3, 4, 5, . . . }. Then J(a, b) = b for all b (see (j^U)) 
and [mm]), i'a,b{x) — {a ~ l)x —: i^a{x) is the same for all &, and denoting Ti = Ta^b and 
T2 = Ta,b' we can replace estimates (|22l) - (|27|) by 



i'^T2 ) "0X2 dm- P'^^ {tpT, hr, ) -0Ti dm 
I J I 



^ J^M^)MT2x)dx- J^M^)MTix)dx (45) 
= (a - If Q'^ X Tl,, (x) dx- jx TZ,{x) dx^ 

To evaluate this difference assume henceforth that < ^ < ^. There is no loss in doing so, 
because Ta^b and Ta,-b are conjugate as obseved above. Define the "rotation" R : I ^ I by 
R{x) = x — mod (Z — i). It conjugates Ta^b to Ta,o, namely 

R{Tl^x) = TJ^^biRx) for all a; e / and n e N. 
Therefore, denoting b = + and Xb{x) = i)^^^ ^ 



xTZbix) dx = / Rix) RiTZox) dx ^ [x + Xb{x)){TZoX + Xh{TZox)) dx 
I J I J I 



so that 



xT:^,,ix)dx- J^xT:\ix)dx 

K.qX ■ iXb'ix) - Xb{x))dx + / PafliXb' - Xb){x) -xdx 
I J I 

+ J^Pa,QXb'{x) ■ {xb'{x) - Xb{x))dx + J^Pa.oiXb' - Xb){x) ■ Xb{x) dx 

As JjXb'{x)dx = and Var(x&) = 2 for aU 6, the third term is of order 0{{2/a)"'\b' - b\) by 
Lemma [3l As P^qX — a^"x, the first term is at most of the same order. Therefore their sums 
over all n are of the order 0(16' — b\). 

We turn to the two remaining terms. Their sum from n = to 00 is of the form 



oc „ oc „ 

/ PaAXb' - Xb){x) ■ gb{x) dx^Y ixb' - Xb)ix) ■ gb{T^flx) dx 



(46) 



with gb{x) = X + Xbix). Let S = ^3^. As J^{^/2(Xb' - Xb){x) dx = 0, Var(x6' - Xb) < 4, and 
\gb\ < 2, the n-th integral is of order C'((2/a)"). Hence the sum from n — Ng [ '"fjg^ 1 to 
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oo is of order \6\, and it remains to estimate the sum from n = to Ns — 1. For these n we 
have a"\5\ < 1. 

We start with the special case 6 = where we have xt = 0, so gb{x) = x. Then the n-th 
term in the sum (j46p evaluates to 

S-{-^ + ^a"6) if a is odd 

^ 2 



5 ■ (ia"(5) if a is even and n > 1 



It follows that the sum from n = to A^a — 1 in (|46)) is of the order Od^l) if a is even and that 
it is (5 • (-^^^ + 0(1)) if a is odd. 

Consider next the case 6=^. As 6 = —^ gives rise to the same map, the we may assume 
w.l.o.g. that -| < 5' < b, i.e. S < 0. If a is odd, then Ta,o6 = b + ^ and Tafi{b + i) = S. 
Therefore 

(Xb' - Xb)ix) ■ gbiK.ox) dx^ 9b{Ta.Qx) dx 

J'b 

n is even 

~ I /o^(^ + 5 + «"^) + ^''(^ + ^ + if " is odd 

_ J (56 + ia"(52 ^ 5 j_A_ if „ is even 

~ \5{b+^) + ^a"S^ -S-^ if n is odd 

as long as a"|(5| < ^^^j, i.e. n < Ng :— Ns — '°^''°iog"'a For the remaining n this identity 
needs to be modified by at most \6\. In any case, 

E / (^'" - ^'') (^) ■9biT-o^)dx = ^-U^-0 + OiS)=S log I ^ TT^ + O (S) for odd a 

with a constant in "O" that depends on 6 and a but not on 6'. If a is even, then Ta,ob — 6, and 
following the argument for odd a and even n we obtain 

Ns-l ^ ^ ^ 

E {Xb'~Xb)ix)-gbiT:^ox)dx = Ns--+OiS)^S\og\S\-\ +0(6) for even a. 

We turn to more general parameters < 6 < i. We have to estimate 



s(^) E / (^''' - ^'')(^) • 5fc(7^a"o^) = E / abiK^.x) dx 

n=0 "'^ n=0 

= E / 9b{Koib)+a"t)dt 

. n "'0 



(47) 



The details of this estimate depend strongly on the distributional properties of the orbit of 6 
under Ta.o and we discuss only two particular but important cases where the situation does not 
become too complicated. First we look at such b for which the orbit of 6 is eventually periodic 
but where the periodic part does neither contain — ^ nor 6. (This is a countable dense set of 
parameters.) In this case one can argue as above for 6 = i and odd a and show that 

s{5) = 5\og\5\-'^ + 0{5) 
log a 

with Ch = Jj gb{x) dfj,b{x) where /if, is the equidistribution on the periodic part of the orbit of 
6. Exceptionally this may be zero, but typically it won't. 
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Next we look at Lebesgue typical points 6, i.e. at Lebesgue typical parameters b. For fixed 
S we interpret s{S) as a random variable where randomness is introduced via the parameter 

b G (0, ^). We are going to show that the random variables S~^Ng ^ s{6) converge in distribution 
to a mixture of Gaussians, i.e. 

CiS-'^NgKiS)):^ [ Af{0,crl)dz as (5^0 (48) 

with suitable variances > that depend on the fixed parameter a. This shows that approx- 
imately s(S) — <5(log Z with a random variable Z that is a mixture of Gaussians which 
depends only on the fixed integer parameter a. 

As a first step we compare s{6) to J2n=o^ ^'9b{Ta.o^)- 9b{^) ~ x + Xb{x), we can estimate 
the difference for the x- and the Xf,(a;)-constributions separately. For the ^-contribution the 
difference is easily seen to be of order 0{6). For the Xb(x)-contribution we estimate each of the 
last Ls l^^^ log A^^] terms of the sum by 2 thus getting a contribution of order 0{S log Ns). 
For the remaining terms we note the following two estimates which are obvious from a short 
look at the graph of T^q: 

m{bel: |T,"o(&) - (-^)| < m{b e I : |r,"o(6) - 6| < a^S} < 4a"<5 , 



so that 



J^bel: 3ne{0,...,Ns-Ls-l} s.th. \T:^,{b) - (-i)| < a^S or |T„"o(6) - b\ < a^S^ 
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- a- 1 a-l ^ 

It follows that Enio"^"'^ Iq Xb{T^^o(J>) + a^t)dt = EnLo"^''^ ^XhiT^ob) except on a set of 
b of Lebesgue measure at most ^^Af^'^. Hence, observing that LgNg'^ as (5 — > 0, the 
convergence in (j48p will follow once we have proved 

C{Yns)^ [ Af{Q,al)dz as Ns ^ oo (49) 
Jo 

where Y/v(6) := A^^^ '}ln=o 9b{T^^b) and 6 is uniformly distributed in the interval (0, i). (To 
ease the notation we abbreviate T^.o by T.) As the single contributions to the sum in Yiq 
depend on b via T"{b) and b itself, this is not the situation of the usual central limit theorem, 
so we treat the problem in two steps: 

Step 1: For fixed z e (0,1) consider Y^{b) A^-2 E„=o 5^ It is a well known general 
fact that, for fixed z, the converge in distribution to some A/'(0, tr^) (see e.g. [2n[T8l[38] ') 
- except for the strict positivity of cr^. To prove this we use ,38, Lemma 6]: suppose for a 
contradiction that cr^ = 0. Then there is a function ijj : I ^ 'E^ oi bounded variation such that 
g^{x) = iPiTx) - ^]{x) for Lebesgue-a.e. x. Let M := 2sup|V'|. Then \ J2l=o gz{T''x)\ < M for 
all n and a. a. x. Looking at suitable periodic orbits it is easy to see that there are n = no and 
X — xq for which this sum is larger than M + 2. But then, as both T and also are at least 
one-sided continuous, there is a small interval close to xq on which the same sum is larger than 
M + 1, which contradicts the above bound that holds for all n and Lebesgue-a.a. x. 

Step 2: We need a number of preparations: 

(i) Let J :— (— i,— i + 2{a-i) ) the interval through which b ranges when b is chosen 
randomly from (0, ^). 

(ii) Let {rj)j > be any sequence of natural numbers tending to infinity and such that rj < 
for all j. For each j denote by Cj C / = [—5,5] a set of points that subdivides / into 
a^^ intervals of the same length which are mapped onto / bijectively by T''^ . (If a is odd 
take Cj := T-''^{-^}, if a is even take Q := T-''^{0}.) For z € Q denote by /| the 
subinterval with left endpoint z. 
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(iii) \\Yn~YnoT'^'''^\\^ = 0{rNN-2) = o(l) as TV ^ oo because the two sums involved differ 
only by 4rjv terms. 

(iv) Var(P^"(;;^lj - 1)) < ilY- y^ri^lj - 1) = o(l) as TV ^ oo by Lemmal 

(v) m{x e /f : Y^iT'-^x) ^ rjv(T'~"a;)} = a-'-"m{a; G / : Y^{x) ^ Yn{x)} < a-'-'^m{lf). 



In order to prove ([38]) we now proceed as follows: it suffices to show that for each bounded 
Lipschitz function : M R holds 

^ ^ (t){YN{b))db^ (^j^(f){b)dN[Q,(jl){b)^dz + o{l) asiV^oo. (50) 



m(J) 

To simplify the notation we write (^O^n) dm instead of Jj (f)(YN(b)) db etc.. 

^ (l){YN)dm:^ [ (f){YN oT^-^^jdm + oil) hy ^ 



J 



m{J) J J m{J) 

P^"{m{J)-Hj)-(l){YNor'^)dm + o{l) 



I 



0(yAror''")dm + o(l) hy ^ 
y / 0(rj^oT''")dm + o(l) by© 

J2 I Pt"'^i^ ■ HYn) dm + o(l) 

J2 MI^) I ^{Y^)dm + o{\) 

= y m(/f ) / cj^dNiQ, ^l) + 0(1) by step 1 

(j)dM{fd,al)] dz + o{l) as N ^ oo 







where one has to choose a sufhciently slowly growing sequence (rjv) in the second last equality. 
Summary of results for integer a In view of the factor (a — 1)^ in ([45]) and of the definition 



of (5 = - — r the above discussion shows: 

a— 1 



b'-b 

(1) For even a > 4, 



Da{b') - Da{0) = 0{b') and 

1^.(6') - 1^4) = 4) log 1^' + - i) 



(2) For odd a > 3, 



Da{b') - DaiO) = -:^r^&' log \b'\-' + Oib') and 
2 log a 

Daib') - Dai-) - ^^—^ib' - i) log 16' - -\-^ + 0(5' - i) 
^ ^2^ 41oga^ 2^ 2' ^ 2^ 
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(3) If b is such that 6 = — ^ + is eventually periodic under Ta,o and the periodic part of 
the orbit neither contains — ^ nor b, then there is a constant Ca.b such that 

Da{b') - Da{b) = C\b {b' - b) log \b' -b\-^+ Oib' - b). 

This generalizes observations (HI and See also Remark O 

(4) For fixed 5 and random b drawn uniformly from (0, ^) or from (— i, 0), 

Da{b + Sia - 1)) - Da{b) = Ca S{log \S\-'y/^ Zs + 0{S) as 5 ^ 

with a constant Ca > and random variables Zg which all have the same distribution - a 
mixture of Gaussians as in (|48p depending only on the fixed parameter a. 

In view of these findings the graphs of Da ■ [—5,5] ^ R are fractal in the sense of Section 12.2 
in [U] - at least for integer values of a -, although they have box - and Hausdorff-dimension 1. 

Remark 4 In Refs. [27l[T4l[32] it has been shown that for 6 = the dynamics of ipx can 
be expressed in terms of generalized Takagi (or de Rham) functions. Analogous conclusions 
hold for the case of & 7^ [S^. The above results are thus intimately related to continuity 
properties of this class of functions under parameter variation. These functions are defined 
by simple functional recursion relations and have been introduced in the literature completely 
independently from the diffusion problem considered here. 

5 Numerical results 

Guided by the analytical results of the previous sections, in this part we numerically study the 
transport coefhcients generated by the piecewise linear maps (I34p . We are aiming particularly 
at a numerical verification of Proposition [2] and of the summary of Section [4.2l To some extent 
the numerics enables us to go beyond the analytical findings as far as detailed local properties 
of these transport coefhcients are concerned. 

Let us start with a reminder of previous results: Exact analytical solutions for drift J and 
diffusion coefficient D for all parameter values a, b were derived in [15] H In [30] data sets were 
generated from these expressions and analyzed by standard numerical box counting |41| . This 
procedure relies on the assumption that 

N{e) ^ e-^ (51) 

for small enough e, where N is the number of square boxes of side length e needed to cover 
the graph of J or D, and B — dims (graph) defines the box (counting) dimension. Analysing 
D{a) = D{a, 0) on 2 < a < 8, see Fig. [1] (a), based on 10^ data points uniformly distributed in 
a yielded a box dimension of S ~ 1.039 [30 . The inset in Fig. [T] (a) depicts N{e) for the new, 
larger data set of 10^ points of D{a) in comparison to (|5T| with the above exponent. Data and 
fit are undistinguishable. 

Fig- [11(b) displays the numerical results from [30] for the local box dimension B{a) of D{a). 
That is, according to ([5T|) B was computed locally on a regular grid of small subintervals Aa 
centered around a. The figure shows that (|5ip yields locally different results for B{a) forming 
an oscillatory structure that becomes more pronounced the smaller Aa. Consequently, D{a) 
was said to be characterized by a "fractal fractal dimension" in [3D] . 

^Another set of formulas was reported in [9] but only for D(a). 
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2 3 4 5 6 7 8 

a 

Figure 1: (a) Diffusion coefficient D{a) on the interval 2 < a < 8 with 2000 data points for 
the one-dimensional map ([34]) sketched in the lower right edge with bias b — and a as the 
slope of the map. The inset of (a) shows standard box counting based on 10^ data points for 
D{a), where N is the number of boxes of side length e. The dashed line in the inset depicts 
the power law ([5T|l with box dimension B — 1.039 as computed in [3^. (b) displays the box 
dimension B{a) computed locally on a regular grid of subintervals of size Aa = 0.006 centered 
around a. For each subinterval a data set of 10® values has been used, and a running average 
was performed over any three neighboring B(a). Both figures, except the inset, are from [30j . 

5.1 Box counting for the diffusion coefficient 

Motivated by Proposition [2] and by [35j , the numerical results of fSU] are now reevaluated 
and supplemented by new, further numerical analysis. We start with the diffusion coefficient 
D{a). Corollary [2] states that B{a) — 1 for all intervals Aa, which is at variance with the 
results presented in Fig.[TJ However, in contrast to the standard box counting assumption (I5ip . 
Proposition [2] is compatible with the existence of multiplicative logarithmic terms by giving 
upper bounds for their exponents. The discussion in Subsection 14.21 shows that these terms do 
indeed exist. 

In detail. Corollary [2] states an upper bound for the box counting function N{e) of D{a) of 

N{e) < Kie-\l-lne)^ . (52) 

This motivates us to plot the product Ne as a function of — In e: For small enough e and in 
double-logarithmic representation one should then see a straight line with the slope yielding 
the exponent of the logarithmic term. Fig. [2] numerically verifies the existence of this term for 
D{a) on 2 < a < 8: There clearly exists a non-zero exponent, however, in the numerics — Ine 
is not large enough to overcome the additive constant in ([52|) for producing a straight line. 

In Fig. [2] three data sets have been plotted consisting of different numbers of data points for 
D{a). The bending off of the graphs at larger — Ine reflects that box counting starts to resolve 
the single points of all the underlying data sets: From the figure one can roughly estimate that 
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Figure 2: N is the number of boxes of length e needed to cover D{a) shown in Fig. [T] (a) 
generated from 10^, 10^ and 10^ data points. Motivated by ([52]) . in contrast to the inset of 
Fig. H] (a) here we plot the product Ne as a function of — Ine double-logarithmically. The 
dashed black line represents a three-parameter fit for the largest data set over 0.5 < — Ine < 12 
with the functional form of ([55)1 . The inset shows results for the exponent of the logarithmic 
correction a obtained from fits where we vary the upper bound emax of the fit interval. 

for a data set of 10^ points for D{a) deviations set in around — In ecut — 7, or ecut — 
Compared with a separation oi da = 6 - 10^^ between any two data points along the a-axis, this 
yields a difference of about three orders of magnitude. The same order of magnitude argument 
holds if one compares ecut obtained approximately for 10® data points from Fig. [2] with the 
corresponding separation of = 6 • 10~* between any two data points. This leads to the 
prediction that for the set of 10^ data points — \necut — 14 in Fig. [2l 

Inspired by (|52|) . we now fit the box counting results with the function 

N{e)^K5e-\l + Ke\ne)°' (53) 

instead of ([51]) . If this fit function reproduces the numerically computed N{e) reasonably well, 
Proposition [2] predicts that < a < 2. However, we emphasize that this Proposition only gives 
us a strict upper bound - it does not actually tell us the "true" functional form of the whole 
graph. We have indeed checked that fit functions others than (|53p . which also obey ([5^ . work 
similarly well. In order to be close to Proposition [2] we stick to the fit function ([53]) in the 
following. 

The dashed black line in Fig. [2]shows a fit of the box counting results for 10^ data points 
of D{a) with this functional formjj The inset of Fig. [2] depicts results for the exponent a 
computed from different fit intervals [0.5, — Inemaa;] for the same data set of 10^ points. It 
indicates convergence towards a ~ 1.2 {—lne,nax 12). The decrease for —Incmax > 12 is 
well in agreement with the cutoff predicted above, which is due to the limited data set. Note 

'''For all fits the nonlinear least-squares Marquardt-Levenberg algorithm as implemented in gnuplot 4.0 has 
been used. 
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Figure 3: Local variation of the product Ne needed to cover D{a) around integer values of a: 
(a) shows results for parameter intervals of size Aa = 0.06 centered around different a, based 
on 10* data points, (b) displays results for subintervals Aa all centered around a — 4, whereas 
in (c) all subintervals converge towards a = 5. In (b) and (c) the graphs have been scaled by 
multiplying e with the order of magnitude difference between the different values for Aa. 



that the cutoff sets in much later than the beginning of the plateau. Hence we conclude that 
for a data set of 10^ points for D{a), 2 < a < 8, and by assuming the fit function (|53p . the 
numerical value for the exponent of the logarithmic term is a ~ 1.2. This is again in agreement 
with Proposition l^PI Note that fits by do not tell the full story: The numerically exact data 
in Fig. [5] show the existence of a non-trivial fine structure pointing towards more complicated 
functional forms for the "true" N(e), which should reflect the intricate structure of D{a) in 
Fig.[l](a). These irregularities may not be understood as numerical errors. 

After having verified the existence of logarithmic contributions on large parameter intervals 
we now look at local variations of the exponent a. This is demonstrated by doing box counting 
for D{a) on small intervals around integer values of a. Fig. [3] (a) reveals that there exist two 
families of curves: The one for even a is at the bottom of this figure, whereas the one for odd 
a is on top. Additionally, all graphs show up such that the ones for larger slopes are always on 
top in both groups thus creating an oscillatory structure. 

We first consider the special case a = 2.03, where according to Fig. [3] (a) a ~ 0. Note that 
D(2) = 0, correspondingly the parameter region just above a = 2 marks the onset of diffusion, 
cf. Fig. [1] (a). As described in [27, 28, 32 , for a ^ 2 there is asymptotic convergence of D{a) 
to the simple random walk solution D{a) — (a — 2)/ (2a). This physical argument explains 
why a — > (a — > 2). There is a trend that larger even integer slopes in (a) give < a < 1 
whereas odd a give 1 < a < 2. Unfortunately, the fits producing these results are very unstable, 
hence even these rough estimates should be taken with care. In any case, the indicated order 

*We have checked that these fit results do not significantly depend on the choice of the initial seeds for 
our three fit parameters and that the asymptotic standard error for them is less than 10% for — In emax > 10. 
However, in our view quantitative error estimates are not reUablo in this case, because we may not assume that 
the residua are normally distributed random variables. 
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of magnitude of a appears to be in agreement with Proposition[21 Our fits furthermore suggest 
that not only a is a function of a but also that the other two parameters in (|53p are locally 
varying. This agrees with conclusions drawn in ^Hj- 

Figs. [3] (b) and (c) provide a more detailed local analysis by looking at successively smaller 
subintervals around two specific slopes. While (c) suggests a ^ (Aa 0) around a = 5, (b) 
with a — A yields approximately a — > 1 (Aa ^0)0 Note that the graphs in (b) and (c) have 
been scaled as described in the figure. Interestingly, this transformation leads to a collapse 
onto a master curve in (b), whereas it does not work that way in (c). Similar observations 
have been reported in [35]. Together with the analytical results of Subsection 14. 2i Fig. [3] thus 
demonstrates remarkable continuity properties of D(a, b) around integer slopes, which strongly 
depend on the direction in parameter space. 

These differences between graphs for odd and even a are consistent with the local box 
counting dimension displayed in Fig. [T] (b): They suggest that the oscillatory structure in 
B{a) actually reflects local variations of the parameters in (|53[) determining the logarithmic 
corrections, erroneously being fit in [30 with the standard box counting equation (|5ip instead 
of taking the existence of logarithmic terms into account. This result is confirmed by covering 
small parameter regions around a — A and a = 5 with non-overlapping sequences of subintervals 
and looking for local variations of the box counting results. Again, one finds oscillations that 
roughly correspond to the ones in Fig.[T](b). Although there is no linear functional relationship 
between a{a) and B{a), one may thus argue that Fig. [1] (b) tells us something about the 
magnitude of local logarithmic corrections. 

We remark that with the computing power available to us it was impossible to produce a 
graph like Fig.[T](b) for local values of a, because for each a{a) box counting would have required 
a data set of at least 10^ values of D{a). Such large data sets appear to be necessary because 
of monotonicity of the exponents: if i? is a subset of F then a{E) < a{F). Local variations 
of a thus pose a serious problem to any numerical box counting analysis, since eventually a(a) 
should always converge to the largest local exponent. However, if this exponent is exhibited 
just on a tiny subinterval it could be extremely tedious to detect it numerically. This argument 
of course also applies to our previous result of a ~ 1.2 for D{a) on 2 < a < 8, which strictly 
speaking only holds for the given data set of 10^ points. We cannot exclude that some tiny 
interval of D{a) eventually yields a larger value of a. In other words, the goal of our numerical 
analysis cannot be to compute unambiguous values for any exponents but rather to demonstrate 
qualitative and quantitative order-of-magnitude agreement with Proposition [2l 

5.2 Box counting for the drift 

We continue our numerical analysis by investigating the parameter dependence of the drift, 
or current, J{a,b). As for the diffusion coefficient, we start with a brief reminder of previous 
results in form of Fig.[4j Like Fig.[Tl it displays a highly oscillatory structure both in the drift as 
well as in the local box dimension as functions of a for fixed 5, where B(a) has been computed 
according to (I5ip by again disregarding any logarithmic corrections. Note particularly the 
pronounced minima at odd integer values. As before, we now reevaluate these findings on the 
basis of Proposition [5] by taking logarithmic terms into account. 

Fig. [5] numerically confirms the existence of logarithmic corrections for J(a,6): There exist 
non-zero exponents a as allowed by Proposition [2l Note particularly the pronounced, different 
fine structures of both curves displayed in the main part, which are much stronger than in Fig. [5] 
for D{a). Due to these oscillations, in case of J{a,b) it is numerically very difRcult to extract 
reliable values for the exponents a by using ([55]) . The two fits included in the main graph yield 
an order of magnitude of a ~ 0.1, which matches to Proposition O 

The inset of Fig. [5] is analogous to Fig. [3] (a) in that it shows box counting results for the 
current J(a, b) , b ~ 0.01, mostly at integer values of the slope a. Note that J(2, 0.01) ~ [15], 
which marks the onset of the drift. As we have argued for the diffusion coefficient, at a = 2.03 
we are thus in a random walk regime for which one may expect a ~ 0, as is shown in the figure. 

^Again, the fit results are highly unstable, so the latter value should be taken with care. 
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Figure 4: (a) Drift J{a) = J(a, 0.01) on the interval 2 < a < 8 based on 2000 data points. As 
in Fig. [U (b) depicts the box dimension B{a) computed locaUy on a regular grid of subintervals 
of size Aa = 0.01 averaged over any three neighboring points. Both figures are from [30] . 

However, a{a) ^ in all the other cases of the inset suggesting again a local variability of a for 
J {a, b), at least around integer values of a. As in Fig. [3] (a) there exist two family of curves, one 
for even a at the bottom and one for odd a on top of the figure. There is also again an additional 
ordering, however, here it is such that curves for larger slopes are always at the bottom in both 
families of graphs, except at a = 2.03. The additional graph for a = 7.09 exemplifies the strong 
local variability of a around a = 7 which, as well as the difference between odd and even slopes 
for box counting results of the drift, agrees with the oscillations in the local box dimension B{a) 
shown in Fig. [4l 

Fits for all the inset curves yield a trend towards small exponents around even and somewhat 
larger values around odd slopes with an order of magnitude of < a < 1, which appears to 
be consistent with Proposition [2] However, we emphasize again that these results give only a 
rough indication for the numerical reasons discussed above. Exact results are only available 
for special cases: As we have discussed in Subsection 14.21 J{a,h) = b for constant a S N under 
variation of b, where we thus have a = 0, linear response and a caricature of Ohm's law. For 
general a one finds that J(a, b)/{b\ log |6||) is bounded but has no limit for 6 ^ [TS] pointing 
towards logarithmic corrections. 

We have also qualitatively checked graphs of D{a, b) and J(a, b) for other parameter values, 
that is, by choosing different values for a and b fixed in the parameter plane and studying the 
resulting functions of the remaining free control parameters. Qualitatively, we obtain results 
that are analogous to the ones discussed above. 

5.3 Continuity properties of the diffusion coefficient at integer slopes 

The previous two subsections demonstrated a very peculiar behaviour of local box counting 
results for drift and diffusion coefhcient around integer slopes a at fixed values of the bias b. 
Subsection 14.21 in turn, gave exact analytical expressions for the difference Da{b') — Da{b) of 
the diffusion coefhcient as a function of A6 = 5' — 6 at integer a in the limit of small A6. This 
suggests to numerically study the continuity properties of Da (b) at fixed integer values of a in 



In order to access suitably small values of the parameter Ab, we have employed the Fortran90 
library mpfun90 pLj for arbitrary-precision arithmetic. Using this library we have calculated the 



more detail. 
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Figure 5: Main graph: Product Ne as a function of — Ine for the drift J(a, h) over the interval 
2<a<8at6 = 0.01 (red '+' symbols, based on 10^ data points) and at 6 = 0.49 (green 'x' 
symbols, based on 10* data points). Included are two fits over the intervals 8 < — Ine < 13 
(6 = 0.01) and 4 < — Ine < 11.5 (6 = 0.49). Inset: Local variation of the product A^e as a 
function of — Ine for parameter intervals of size Aa = 0.06, mostly centered around integer 
values of a and based on 10® data points. The graph for a = 7.09 demonstrates that, in 
agreement with Fig. [4] (b), there exist strong local fluctuations of the box counting functions 
under variation of the slope a of the map. 



difference quotient {Da{h') - Da{b))/{h' -b) oiD with fixed 6 at values of A6 down to 10~200^ 
Figure [6] (a) shows a subset of our results for a = 3 and a = 4 at fixed b e {—0.5, 0}. There is 
excellent agreement between the numerical results and the analytical observations ([1]) and ^ 
of Subsection 14. 21 predicting straight lines. This agreement is as good to the limits of attainable 
precision, and has been checked for other integer values than those shown in Fig. [6l 

Figure[n](b) depicts the diffusion coefficient Da{b) at a = 4 and a blowup around 6 = 0, which 
corresponds to the two curves in (a) at this a value. Note that there is reflection symmetry 
for Da{b) with respect to & = —0.5 and b — 0. One can see that at & = —0.5, where the 
difference quotient in (a) displays a multiplicative logarithmic term, Da{b) in (b) exhibits a 
global maximum in form of a sharp cusp. The global minimum at 6 = 0, on the other hand, is 
approached in a rather smooth, oscillatory manner yielding a rounded-off shape, see the inset in 
(b). This relates to the difference quotient curve in (a) with zero logarithmic term. Analogous 
observations are made for a = 3, where Da{b) exhibits local maxima both at b = —0.5 and at 
6=0, and for other integer slopes. We remark that the quite regular structure of Da{b) in (b), 
particularly around both local extrema, resembles very much the one of the fractal generalized 
Takagi functions studied in ^27|.I14|[32]. 

Observation ^ generalizes observations ([T]) and ^ by stating that logarithmic corrections 
are typical for parameter values of b yielding Markov partitions. In [26l|27l[29l[T5] it has been 
shown (for 6 = 0) that Markov partition parameter values identify local maxima and minima of 
the parameter-dependent diffusion coefficient by relating them to ballistic and localized orbits 
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Figure 6: (a) Difference quotient {Da{b') — D{h))/{h' — 6) as a function of ln(6' — h) at integer 
values of a e {3,4} and with fixed h G {—0.5,0}. Each curve is based on 10^ data points. 
Included are best fit curves (black dashed lines) whose fitted slopes (from bottom to top: 
3/(2 In 4), l/ln3 and 1/(2 In 3)) agree with the analytic predictions of Subsection 14.21 to four 
significant figures. The case a = 4, 6 = clearly has slope zero, as predicted. The barely visible 
fine scale oscillations of each curve reflect higher order correlations in these quantities, (b) 
Diffusion coefficient Da{h) at a = 4 for —0.5 < 6 < and a magnification of the region around 
6 = 0. For each curve 2000 data points have been computed from exact analytical solutions for 
Da{b) [T5^. These curves form the basis for the two graphs at a = 4 displayed in (a). 



of the critical points of the lifted map, respectively. One may thus speculate that the above 
numerical observation holds true for local extrema on finer scales, that is, that local cusps in 
Da{h) reflect logarithmic corrections in the local difference quotient, whereas rounded-off local 
extrema signify the lack of logarithmic terms. See also [35] for related results. Furthermore, in 
Fig. [6] (a) we deliberately restricted the range of ln(&' — h) so that, upon very close scrutiny, a 
fine structure of all curves can be seen on top of the straight hue behaviour. Fig. [6] (b) suggests 
that this oscillatory fine structure, which yields higher order corrections to the analytical results 
of Subsection 14.21 is induced by the fine structure of Da{h). 

We have also numerically investigated the accuracy of observation ^ in Subsection 14.21 
Its main statement is that at fixed A& and with h values taken uniformly from the interval 
[0, 1/2), the quantity {D{b')~D{b)) / Ab\/— In A6 should be distributed like a mixture of centered 
Gaussians, that distribution being independent of the particular value of Ab. In fact what is 
typically seen at integer slopes is a distribution rather close to a pure Gaussian. We have tested 
this using the technique of quantile-quantile plotting (qqplots) as well the standard Shapiro- 
Wilk normality test. Both tools were implemented in the statistical package R [42 . 

Figure [7] presents results obtained for three sets of data with the slope fixed at a = 4. For 
larger a, the results become closer to a fixed Gaussian, as the function g{x) in (|46|) becomes 
more dominated by the x term which has no b dependence. Here however, deviations from 
Gaussianity can be seen, at least for sufficiently small Ab. In the three parts of Figure [3 the 
red line with slope a and zero offset ^ shows the theoretical result for a Gaussian distribution 
with standard deviation a and mean fi, with those parameters here taken as those of our data 
set. As can be seen, all our distributions show close agreement with this curve. However, the 
Shapiro- Wilk normality test is more discerning: in (a) Ab — 10^^° and we obtain a p- value of 
only 0.008, well below the significance level for rejecting the null hypothesis of normality. In 
(b) Ab = 10"^" and we get a p-value of 0.25, demonstrating that this distribution is indeed 
very close to a pure Gaussian. It is however likely that the deviations from Gaussianity in (a) 
are rather due to deterministic effects arising from the relatively large value of Ab chosen and 
not from the nature of the true limiting distribution being a mixture of Gaussians predicted 
by observation ([4]) . It seems that the dominant behaviour when the distribution seems to have 
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Figure 7: (a) Normal quantile-quantile plot at a = 4 for the distribution of D{b') — D{b), scaled 
by A&V— In A6, with chosen Ab = b' — b = 10~^" held constant and b picked from a uniform 
distribution on [0, 1/2). The red line with slope a and zero offset ^ would be the result for a 
Gaussian distribution with standard deviation a and mean /i. Here the numerically obtained 
values of these parameters were used for the fit. (b) As (a) but with A6 = 10~^°. (c) As (b) 
but with the range of b restricted to [0,0.005). 

converged is not detectably different from a pure Gaussian. We note that despite this, the two 
distributions in (a) and (b) are similar and both have mean close to zero, demonstrating that 
we have no disagreement with observation ([4]) , merely that its details are too sensitive to check 
numerically. 

It is however possible to go further numerically, for example one can also study the nature 
of the distribution obtained when b is taken from a subinterval of [0, 1/2), which as can be 
seen in Fig. [7] (c) leads in the case of integer a to distributions with rather more fine structure 
than the nice curves seen in Figs. [7] (a) and (b). This is clear evidence of the deterministic 
nature of the underlying system in the form of strong correlations at fine scales. In this case 
the Shapiro- Wilk p-value is about 0.001. 

As far as observation ^ is concerned, away from integer values of a quite different behaviour 
is seen thus clarifying that this observation is rather to be considered atypical. Here the 
distribution of 13— differences is centred around zero still, but with a more sharply peaked and 
heavily tailed distribution than a Gaussian. These deviations persist even very close to the 
integer cases (e.g. at a 3 + 10^^°), though Gaussian behaviour does appear to be approached 
slowly in the limit of integer values. 

These numerical methods can also be used to investigate variation of the continuity of the 
transport coefficients as b is held fixed and a varies, as considered in [35] and already looked 
at using box counting in Fig. [3] Here the maximal exponent of logarithmic correction, i.e. 
D(a') — D{a) ^ \a' — a\{h\\a' — a|)^, can be seen for odd a, and though this might appear to 
be in contradiction to the third part of Fig. [3] for a = 5, in fact arbitrarily close to a = 5 the 
exponent tends locally to zero. Thus the box counting only sees the "typical" local behaviour 
and the current method is more suited for picking out atypical behaviour at specific points. 

6 Conclusions and outlook 

(1) We proved rigorously that the diffusion coefficient of deterministic random walks generated 
by piecewise expanding interval maps depends continuously on the maps. More precisely, 
for "natural" parametrizations of the maps by some parameter A, the diffusion coefficient 
as a function of the parameters has a modulus of continuity not worse than |(5A|(log |<^A|)^. 
Even if all maps in the family are topologically conjugate, the detailed analysis of section l4.2l 
shows that the modulus of continuity cannot be expected to be better than \5\ ■ \og\5\\\. 
This is in sharp contrast to the situation for the drift (or other averages of observables) 
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that depend difFerentiably on parameters in this case [1]. One might thus conjecture: If 
the maps are all topologically conjugate as in the case of integer slopes, then J is Lipschitz 
and D has simple logarithmic corrections. Otherwise J has simple logarithmic corrections 
and D has quadratic ones. 

(2) We verified numerically the existence of logarithmic corrections in the box counting data 
for both the parameter dependent drift and diffusion coefficients. The computed values 
for the exponents of these logarithmic terms are compatible with the bounds predicted by 
our mathematical theory. However, we emphasize again the serious difficulties to obtain 
quantitatively reliable numerical results, which required to analyze huge data sets. These 
difficulties are due to strong local variations of these exponents and of the other control 
parameters governing the logarithmic corrections, as we find numerically. 

These new numerical results correct and amend the previous box counting analysis of 
Klages and Klaufi [30' along the lines conjectured by Koza [35]. Our model thus gener- 
ates interesting examples of fractals for which the definition of the standard box counting 
dimension is misleading. We conclude that the (local) non-integer variations of the box 
dimension reported in f30^ actually reflect non-trivial local variations of the parameters of 
the logarithmic corrections. 

We have furthermore numerically verified analytical predictions for the difference quotient 
of the diffusion coefficient as a function of the bias at integer slopes. These results suggest 
that the existence of logarithmic corrections is intimately related to the shape of the extrema 
in the diffusion coefficient curves. 

(3) In [33 a nonlinear generalization of our present model has been studied, which exhibits 
anomalous diffusion generated by marginal fixed points. Computer simulations led to con- 
jecture that the anomalous diffusion coefficient of this map is discontinuous on a dense set 
of parameter values. It would be interesting to check this conjecture mathematically. 

These fractal transport coefficients also seem to provide a nice testing ground for methods 
of multifractal analysis jTT]. 

Another important problem is to check whether such logarithmic corrections in transport 
coefficients might also be expected to occur in more 'physical' systems, which are perhaps 
even accessible experimentally. This seems to be strongly related to the question whether 
a family of physical dynamical systems shares the same topological conjugacy class under 
parameter variation. 
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